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Abstract 

We consider the design of asymmetric multiple description lattice quantizers that cover the entire 
spectrum of the distortion profile, ranging from symmetric or balanced to successively refinable. 
We present a solution to a labeling problem, which is an important part of the construction, along 
with a general design procedure. This procedure is illustrated using the Z 2 lattice. The asymptotic 
performance of the quantizer is analyzed in the high-rate case. We also evaluate its rate- distortion 
performance and compare it to known information theoretic bounds. 

I Introduction 

A multiple description source encoder generates a set of binary streams or descriptions 
of a source sequence, each with its own rate constraint. The transmission medium may 
deliver some or all of the descriptions to the decoder. The objective is to minimize the 
distortion between the source sequence and the decoded sequence when all the descriptions 
are available, while ensuring that the distortion which results when only a subset of the 
descriptions are available remains below a pre-specified value that depends on the subset. If 
there are D descriptions, the distortion profile is a vector of length 2 D whose components 
give the distortion constraints for each subset of the descriptions. 

In recent years, multiple description coders have received considerable attention, driven 
by the interest in packet voice and video communications (see the bibliography). Most of 
the work (with the exception of has centered around the successively refinable case and 
the balanced/symmetric case, which are in a sense two extremes of the distortion profile. 
Successive refinement coders find application in networks with a priority structure whereas 
balanced codes are useful in networks that do not have such a structure, the best example 
at the present time being the Internet. 
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In this paper we propose a structured scheme that bridges the two cases, in the sense 
that it permits a fairly general distortion profile to be specified. By allowing the individual 
descriptions to have different distortions, the quantizer behavior can range from the balanced 
case (where each description is equally important) to a strict hierarchy (where the loss of 
some descriptions could make decoding impossible). The new design is described in terms of a 
lattice vector quantizer, but the general principle of asymmetric multiple description coding 
can be extended to many other quantizers, such as trellis coded quantizers, unstructured 
vector quantizers, etc. This could potentially allow us to incorporate channel (or network 
route) reliability information into the transmission. Also, it might be a useful way to allow 
for less intrinsic wastage of network traffic as some descriptions could be given to the decoder 
without necessarily waiting for the more important descriptions to arrive (as in successive 
refinement). 

For previous work on the information theoretic aspects of the multiple description prob- 
lem see H [10, |22], |29|, [30]]. The problem of designing quantizers for the multiple description 



problem has been considered in [JTTJ, [17], [19], ^5], |26|, ^8|. The work presented here extends 
that in | 8], which considered only the balanced/symmetric case. Unlike the work in |p.l]| . 



we do not use a training approach; instead we use the geometry of the underlying lattice to 
solve a labeling problem. Other approaches to multiple description coding based on over- 
complete expansions are presented in |1], ^ [L3[] and methods based on optimizing transforms 
and predictors are presented in |16], [21], p7[ . 

The paper is organized as follows. The source coding problem is formulated in Section 
[TI], the design method is described in Section [Til], properties of the lattices and sublattices 
needed for the construction are developed in Section [TV], a high rate analysis is presented in 
Section M, and numerical results are presented in Section [VI. 



II Preliminaries 



Channel 1 



Source 



Q 



Channel 2 



Decoder 



Decoder 



Decoder 



a; 



Figure 1: Block diagram of a multiple description vector quantizer. 

A block diagram of a two-channel multiple description vector quantizer (MDVQ) using a 
lattice codebook is shown in Fig. 0]. An L-dimensional source vector x is first encoded as the 
closest vector A in a lattice A C M L . We will write A = Q{x). Information about the selected 
code vector A is then sent across the two channels, subject to rate constraints imposed by 
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the individual channels. This is done through a labeling function a. At the decoder, if only 
channel 1 works, the received information is used to select a vector X[ from the channel 1 
codebook. If only channel 2 works, the information received over channel 2 is used to select 
a code vector A' 2 from the channel 2 codebook. If both channels work, it is assumed that 
enough information is available to recover A. 

We will assume that the channel 1 and channel 2 codebooks, denoted by Ai and A 2 
respectively, are sublatticesf] of A. The index [A : Aj] is denoted by Ni, % = 1,2. iVj is also 
called the re-use index of sublattice Aj. We assume that each Aj is geometrically similar to 
A, i.e. that Aj can be obtained from A by applying a similarity (a rotation, change of scale 
and possibly a reflection). To simplify the analysis we will usually assume that the Aj are 
strictly similar to A, i.e. that reflections are not used. 

Property 1 Let A be an L-dimensional lattice with generator matrix G ( the rows of G span 
A). A sublattice A\ C A is geometrically strictly similar to A if and only if the following 
condition holds: there is an invertible L x L matrix U\ with integer entries, a scalar c\, and 
an orthogonal L x L matrix K\ with determinant 1 such that a generator matrix for A\ can 
be written as 



d = U X G = aGKi . (1) 
If (fl|j holds then the index of Ai in A is equal to 



rA . . /detAi detd 
* = ' = V detX = ~detG~ 

= det U x = c{ . (2) 

Furthermore, A\ has Gram matrix 

A ± = dGf = UtGG tr Ul r = U x AUf = c\A , (3) 

where A = GG tr is a Gram matrix for A. 

Even if the similarity is not strict, equations (]T|), (|2|) and (|3]) still hold but with det K\ = 

-1. 

We will also usually assume that the sublattices Ai and A 2 are clean , that is, no point 
of A lies on the boundary of the Voronoi cells of A x or A 2 . Our algorithm still applies if this 
condition is not satisfied, but the book-keeping becomes more complicated. 

Finally, we require a sublattice A s of Ai f] A2 which is geometrically strictly similar to 
A and has index N s = N1N2 in A. To reduce the complexity of the design we will also 
sometimes make use of a sublattice A; cm of Ai H ^2 which has index Ni cm = lcm(Ai, N2) in 



A (such sublattices do not always exist - see Section |V|) . 

Since the information sent over channel 1 is used to identify a code vector Ai £ Ai, and 
the information over channel 2 is used to identify a code vector A2 £ A2, we will assume 
that the labeling function a is a mapping from A into Ai x A2 and that (Ai,A2) = a(A). 

1 Strictly speaking, the codebooks are finite subsets of the sublattices Ai and A 2 , but we will ignore that 
distinction in this paper. 
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The component mappings are Ai = «i(A) and A2 = 02(A). In order to recover A when both 
channels work, it is necessary that a be one-to-one. 

Given A, A 1; A 2 and a, there are three distortions and two rates associated with the 
quantizer. For a given source vector x mapped to the triple (A, Ai,A 2 ), the two-channel 
distortion d is given by ||x — A|| 2 , the side distortions di by ||a; — A.j|| 2 , i = 1,2, where 
|| x ||2 jl [1/L)^2f =1 xf is the dimension-normalized Euclidean norm. The corresponding 
average distortions are denoted by do, d\ and d 2 . (We will also refer to do as the central 
distortion.) We assume that an entropy coder is used to transmit the labeled vectors at a 
rate arbitrarily close to the entropy, i.e., Ri = H(ai(Q{X.)))/ L, i = 1,2, where H is the 
binary entropy function. The problem is to design the labeling function a so as to minimize 
d subject to di < D x and d 2 < D 2 , for given rates (i?i, R 2 ) and distortions Di and D 2 . 

We will assume that the source is memoryless with probability density function (pdf) 
p. The L-fold pdf will be denoted by pl where Pl((xi, x 2 , . . . ,£l)) = Yli=iP( x i)- The 
differential entropies satisfy the relation h{pL) = Lh{p). 

Given a lattice A, a sublattice A' and a point A' G A', we denote by Va : a'(A') the set of 
all points in A that are closer to A' than to any other point in A'. This set is the discrete 
Voronoi set of A' in A. If A' is a clean sublattice of A we do not need to worry about ties 
when calculating Va. : a'(A'). The Voronoi cell Va(A) of a point A € A is the set of all points 
in R L that are at least as close to A as to any other point of A. Also S(X') = a(V\ : \/(X')), 
A' G A', will denote the set of all labels of the points in Va : a'(^0- 



A Distortion Computation 

The average two-channel distortion do is given by 

do = / ||a; — \\\ 2 pL{x)dx. (4) 

Since the codebook of the quantizer is a lattice, all the Voronoi sets in the above summation 
are congruent. Furthermore, upon assuming that each Voronoi cell is small and letting v 
denote the L-dimensional volume of a Voronoi cell, we obtain the two-channel distortion 

do = jv ^ fdX = G(A)u 2/L , (5) 

where the normalized second moment G(A) is defined by (f|) 

/ llxll ^dx 

G(A) = Wi (6) 

When only description i is available, for i = 1,2, the distortion is given by 

^ = rfo + ^ll A -«*( A )f P ( A )> (7) 

AGA 



where -P(A) is the probability of lattice point A, and we have assumed that A is the centroid 
of its Voronoi cell. This is true for the uniform density. For nonuniform densities, there is an 
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error term which goes to zero with the size of the Voronoi cell. The first term in (|7|) is the 
two-channel distortion and the second term is the excess distortion which is incurred when 
only description i is available. Note that, for a given A, only the excess distortion term is 
affected by the labeling a. 

At this point we impose a constraint on the labeling function that allows us to reduce the 
problem to that of labeling a finite number of points. We assume that the labeling function 
has the property that a(X + \ s ) = a(\) + \ s , for all X s G A s . This leads to the following 
simplification: 

d i = d + {l/N a ) ll A -«*( A )f> (8) 

AGVA:A s (0) 

where we have assumed that -P(A) is approximately constant over a Voronoi cell of the 
sublattice A s , but may vary from one Voronoi cell to another. 

B Rate Computation 

Let Ro bits/sample be the rate required to address the two-channel codebook for a single 
channel system[[ We first derive an expression for R and then determine the rates R\ 
and R 2 . We use the fact that each quantizer bin has identical volume v and that Pl{x) is 
approximately piecewise constant over each Voronoi cell of Ai and A 2 . This assumption is 
valid in the limit as the Voronoi cell become small and is standard in asymptotic quantization 
theory. 

The rate Ro = H(Q(X)) is given by 

ft = -(1/L)W M^logJ 

x Jv A (\) JVa(X) 

~ -(!/£) 5^ / p L (x)dx\og 2 p L (X)u 

« h{p)-(l/L)log 2 (u). (9) 

It can be shown that the rate for description % is given by 

R i = R -(l/L)log 2 (N l ), z = l,2. (10) 

A single channel system would have used Rq bits/sample. Instead a multiple description 
system uses a total of R± + R 2 = 2R — (1/L) \og 2 (NiN 2 ) bits/sample, and so the rate 
overhead is R - (1/L) log^N^). 



Ill Construction of the Labeling Function 

Suppose A is an L-dimensional lattice with a pair of geometrically strictly similar, clean 
sublattices Ai and A 2 , and let A s (the product sublattice) be a geometrically strictly similar, 

2 This quantity is useful for evaluating the two-channel distortion as well as for evaluating the rate overhead 
associated with the multiple description system. 
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clean sub lattice of both Ai and A2, with indices [A : Ai] = N\, [A : A2] = Af 2 and [A : A s ] = 

In order to construct a labeling function we first identify £ , the subset of points of 
Ai x A 2 that will be used to label the points of A. Next, a one-to-one correspondence will 
be established between Vaa s (0) and a proper subset of £ so as to minimize an appropriate 
objective function, while ensuring that the labeling can be extended uniquely to the entire 
lattice. To this end we first start by formulating a cost criterion that will be used in the 
design. 

A Cost criterion 

The multiple descriptions problem may be formulated || as a problem of minimizing the 
central distortion subject to constraints on the side distortion. The associated Lagrangian 
cost criterion is given by 

2 

J = do + J^Udi (11) 

2 

= (71 + 72 + l)do + ^Y II A - a ^ H 2p ( A ) 

i=i aga 

2 

= (71 + 72 + iH + ]Tp(A) £ 7< ||a - ai(\)\\ 2 . 

AgA i=l 

where 71, 72 are Lagrange multipliers. 

The central distortion do is determined by the lattice A. If we assume that P(X) is 
approximately constant over the Voronoi cell of A s , we can rewrite the cost criterion in 
terms of the cost over a Voronoi cell of A s . Then the design problem reduces to finding a 
labeling scheme a (A) which minimizes 

±r Yl M A -«i( A )H 2 + 72||A-a2(A)|| 2 ]. (12) 

S AgVa ; a s (0) 

After some algebra, the expression inside the summation can be rewritten as 

7l72 H m mM2 . / , \ 11 \ 7l"l(A) +72«2(A) || 2 . s 

l |a 2 (A) -ai(A)|| +(7i+72)||A ■ \\ . (13) 



7i + 72 7i + 72 

The values of 71 and 72 determine the relative values of the two side-distortions d\ and g?2- 
Therefore our design principle is (informally) for a given pair 7x and 72, to find a labeling 
function a(A) such that the sublattice points «i(A) G Ai, 012(A) 6 A 2 are not very far apart 
and the lattice point A G A that is being labeled is not very far from the weighted mean 
(the second term of ([IB])) of these two sublattice points. This general guiding principle leads 
to our lattice design. We will first describe the basic quantizer design and then illustrate it 
using the lattice Z 2 . 
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B Lattice Quantizer 

The quantizer construction is based on the following steps. 



1. We are given an L-dimensional lattice A, rates Ri,R 2 and distortions Di,D 2 . These 
determine the indices N\, N 2 using (|i"0l), and we attempt to find (strictly similar, clean) 



sublattices Ai,A 2 with these indices, together with a product sublattice A s . We also 
choose appropriate values for the weights 71 and 72. For example, a successively 
refineable quantizer corresponds to choosing 71 = 1 and 72 = 0. For the balanced case 
we take 71 = 7 2 . By appropriately choosing A r 1 , N 2 , 71, 72, one can achieve different 
levels of asymmetry in rate and distortion. 

2. We find the discrete Voronoi set^ V = Va : a s (0) for the sublattice A s . This is the 
fundamental set of points that we will label. The labeling is then extended to the full 
lattice using the shift invariance property (see Section ||). We also find the sets 

Vi = VWO) = vW0)nA a , (14) 
v 2 = v M . A M = K A:As (o)nA 2 , (15) 

which are the points of Ai and A 2 belonging to the discrete Voronoi set. 

3. We determine the set 



A(Ai) = {A 2 e A 2 : A 2 G Vo + Xx} (16) 

for all Ai G V\. These are the points in the sublattice A 2 which are in the Voronoi set 
Vo of A s when translated to be centered at Ai G V±. This ensures that the edge length 
|| Q; 2(A) — ai(A)|| 2 will be minimized (see Property 3). We will show that each member 
of £i(Ai) lies in a different coset with respect to the sublattice shifts in A s (Property 
2). Similarly, we determine the set 

£ 2 (A 2 ) = {Ax G Ax : Ai G Vo + A 2 } (17) 

for all A2 G ?2- The set of edges emanating from Vo is given by 

Sedges = {(Ai, A 2 ) : X 1 G V l , A 2 g C 1 (X 1 )} (J {(A 2 , A x ) : A 2 G V 2 , A x g £ 2 (A 2 )} . 

(18) 

We find a set of coset representatives £q for the equivalence classes of £ edges modulo 
A s . Property 6 will establish that we can write £q either as 

£0 = {(Ai, A 2 ) :XiEVi and A 2 G A(Ai)} (19) 

or equally well as 

£0 = {X 2 G V 2 and Ai G £ 2 (A 2 )} . (20) 



3 We usually omit the word "discrete" when referring to this set. 
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4. Matching the edges to the lattice points in the Voronoi set is now a straightforward 
and easily solved assignment problem (cf. fl8fl). To formulate this assignment problem 
we compute the cost given by ([12|) for each lattice point and each equivalence class of 
edges modulo A s (taking the minimum over the edge class). We use only one member 
from each edge class modulo A s in order for the shift invariance property to be satisfied. 
This allows us to construct the set of edges which will later be used to label the points 
in V A:As . 

If there exists a sublattice Ai cm (as defined in Section II) which is also a geometrically 
strictly similar, clean sublattice of Ai and A 2 the computational complexity of the 
design can be further reduced. For then we need only label the coset representatives 
for 8 modulo A icm . We will show that this does not reduce the performance of the 
quantizer — see Property ||. In this case we replace the sets V\ and P 2 by the sets 
V x = VA i; A, cm (0) and V 2 = VA 2: A icm (0). The rest of the procedure is unchanged. 



C Properties of the quantizer 

In this section we state some of the properties of the construction proposed in Section |B[ 
We have imposed the following restrictions on the labeling scheme: 

1. The labels are invariant under shifts by the product sublattice A s . 

2. The labels arise from different cosets of the product sublattice: i.e. if (Ai,A2) and 
(Ai,A 2 ) are valid edges then A 2 and A 2 are in different cosets with respect to the 
product sublattice. 

Property 2 Each member of Ci(Xi) lies in a different coset with respect to the sublattice 
shifts in A s , and |£i(Ai)| = N±. Similarly each member o/£ 2 (A 2 ) lies in a different coset 
with respect to the sublattice shifts in A s , and |£ 2 (A 2 )| = N 2 . 

Proof: Let A 2 , A' 2 £ £i(Ai), and A' 2 = A 2 + A s for some X s £ A s . Then A 2 — A x = A 2 — Ax + A^. 
Hence A 2 — Ai and A 2 — Ai cannot both lie in Vq. But since A 2 , A 2 £ £i(Ai), A 2 — Ai and A 2 — Ai 
are in Vo, a contradiction. Thus each A 2 £ £i(Ai) is in a different coset with respect to the 
sublattice shifts in A s . Now {Vo + A 1 } A ' =Ai+As x s eA s * s a partitioning of the points of A, and 
each of these disjoint sets contains points from different cosets of A 2 (with respect to shifts 
in A s ). Since there are only N\ different cosets of A 2 , |£i(Ai)| < N\. In fact equality must 
hold, because the space is tiled by such sets and if there were a A x for which |£ i (A 1 )| < N% 
then would be a Ai for which | jCi ( Ai ) | > A r 1 , which impossible. An identical proof holds for 
£ 2 (A 2 ). ■ 

Property 3 £i(Ai) consists of the Ni points A 2 £ A 2 closest to Ai subject to the constraint 
that each A 2 is in a different coset. 

Proof: We know that A 2 £ £i(A x ) A 2 £ V + A x <^> A 2 - A x £ V . Then A 2 - A x £ V <^> 
||A 2 — Ai|| 2 < ||A 2 — Ai + A 2 || 2 , for all X s £ A s . Thus for any A 2 = A 2 + A s , A s ^ we have 
||A 2 — Ai|| 2 < ||A 2 — Ai|| 2 , and the claim follows. ■ 
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Property 4 A 2 G £i(Ai) <^> Ai G £ 2 (A 2 ). 

Proof: For clean lattices, if x G V then —x G Vq f|. Then A 2 G £i(Ai) A 2 G Vb + Ai 4» 
A 2 - Ai G 1/ ^ Ai - A 2 G <£> Ai G V + A 2 ^ A x G £ 2 (A 2 ). ■ 

Property 5 As lattice points in A are labeled, the number of times each point from Ai is 
used is Ni and the number of times each point from A 2 is used is N 2 . 

Proof: Let A r (Ai) denote the number of lattice points labeled by Ai G Ai. Certainly 
N(Xi) > Ni since Ai is used N± times when we form the edges {(Ai, A 2 ) : A 2 G £i(Ai)}. If Ai 
is used in more than A r 1 labels then there is a valid edge (Ai, A 2 ) with A 2 ^ £i(Ai). But this 
is impossible by Property [|. Therefore N(Xi) = N± for all Ai G Ai and similarly N(X 2 ) = N 2 
for all A 2 G A 2 . ■ 

Property 6 The number of cosets in the edge set So modulo A s is equal to the number of 
lattice points in Vq. 

Proof: Consider the edge set £ = {(Ai, A 2 ) : Ai G V\, A 2 G £i(Ai)}. From Property || each 
A 2 G £i(Ai) lies in a different coset modulo A s and hence each edge (Ai, A 2 ) G lies in a 
different coset. As |£o^| = A r iA r 2 , there are at least that many cosets in the edge set. ■ 

Property 7 The labeling scheme produces a unique label for each lattice point. 

Proof: This is immediate from the fact that the labels for the cosets of A/A s are taken 
from distinct cosets of S /A s . ■ 



Property 8 The labeling scheme minimizes the cost criterion given in (|TT|) subject to the 
coset restriction. 

Proof: This is an immediate consequence of Property |3|. ■ 

Property 9 Suppose N% and N 2 are not relatively prime, and there exists a sublattice A; cm 
with index /cm{A r 1 , N 2 } in A which is a geometrically strictly similar, clean sublattice of Ai 
and A 2 , and contains A s . Then we may construct the labeling to be invariant under shifts by 
Ai cm , and obtain the same edge set as if we used the product lattice A s . With this procedure 
it is necessary to label only lcm{ N%, N 2 } lattice points rather than N\N 2 points. 

Proof: If such a A icm exists then we just need to show that the edge set constructed by 
using the algorithm with A s can be produced by sublattice shifts of the edge set constructed 
using A [cm . As we saw in the proof of Property ||, the coset representatives for the edge 
set are constructed by using S = {(Ai,A 2 ) : Ai G Vi,X 2 G £i(Ai)}. However, V\ = 
UwA icm UA ie ?v( A i + Ajcm). Therefore £ = \Jx lcm eA lcm 4 where £' = {(A 1; A 2 ) : X 1 G 
V 1 , X 2 G £i(Ai)}. It follows that there are exactly lcm{ Ni,N 2 } coset leader edges in £ 
with respect to the sublattice A/ cm and they are given in S . Therefore, by matching the 
cosets of the edges modulo Ai cm with the lattice points in the Voronoi set for A; cm , using the 
assignment algorithm (as before), and then shifting by A/ cm we produce exactly the same 
labeling as we obtained using A s . ■ 
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Property 10 If there exist several labeling schemes achieving the same cost we can mix 
these configurations to achieve different levels of asymmetry. A sufficient condition for this 
to occur is for the number of unique representation points to be smaller than the number of 
lattice points in the product lattice. 

Proof: The number of representation points is equal to the number of lattice points in 
the Voronoi set Vq (see Property |6|). Therefore, if there are some representation points 
which overlap (i.e. the number of unique representation points is less than the number of 
points in Vo, then there is more than one labeling scheme that produces the same Lagrangian 
71^1+72^2, with each labeling producing different d\, d 2 . Suppose one extremal configuration 
produces the lowest d™ tn and therefore produces the largest d™ ax and if the other extremal 
configuration produces, the highest d™ ax and the lowest d 2 nm . Then by using the first con- 
figuration in proportion a and the second in proportion a = 1 — a one can produce side 
distortions d\ = ad™ m + ad 7 p ax and d 2 = ad 2 nax + ad! 2 am . Thus by keeping the Lagrangian 
cost the same, one can obtain different levels of asymmetry in the distortions d±, d 2 . ■ 

D Example 

In this section we illustrate the design procedure with an example in two dimensions using 
the lattice Z 2 . We choose |Ai| = 5 and |A 2 | = 9. Portions of the two sublattices are shown 
in Figure ^ where the points of Ai are marked with circles, the points of A 2 with crosses, 
and the points of A s with both circles and crosses.^ There are 45 points in the Voronoi 
set Vo for A s . The set V\ contains 9 points of Ai and the set V 2 contains 5 points of A 2 . 
The edges £ e d g es (see Eq. (|i~8D) emanating from the points of V are also shown. These 
are found using the sets £i(Ai) and £ 2 (A 2 ) for A 1 £ A 1; A 2 £ A 2 . For example, if we 
take the point Ai = (2, 1) £ Vi, we see that there are 5 points in the set £i(Ai), namely 
{(0, 0), (0, 3), (3, 3), (6, 0), (3, 0)}. Note that there are several edges emanating from V which 
are a sublattice A s shift apart. For example the edge {(—2, —1), (—6,0)} is a sublattice A s 
shift away from the edge {(4,2), (0,3)}. To satisfy the shift invariance constraint, we must 
use only one of these edges to label a point in Vq. This constraint is built into the optimization 
procedure. The result of the optimization procedure is illustrated in Figure |3|. Here we have 
shown only the points in Ao- The points in Ai flAo are marked by circles and those in A 2 H Ao 
by crosses. Each point carries a pair of labels (Ai, A 2 ) with Ai £ Ai, A 2 £ A 2 . In this example 
we have set 71 = 9 and 7 2 = 5, which determines the respective distortions di obtained by 
the design. A comparison of these distortions with that predicted by information theory is 
given in Section [VT]. 

IV Good lattices 

The lattices that we will investigate and apply in this paper are Z n for n — 1, 2 or a multiple 
of 4, together with the root lattices _D 4 and E 8 ||. The analysis could be extended to treat 
other lattices such as Z 3 , Z 6 , the 12-dimensional Coxeter-Todd lattice, the 16- dimensional 

4 In the enhanced (pdf) version of this document the circles are blue and the crosses are red. 
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Figure 2: Edges emanating from the Voronoi cell of the sublattice. 



Barnes- Wall lattice or the 24-dimensional Leech lattice 
these here. 



but we shall not discuss 



A The construction of similar sublattices 

We begin with the observation that multiplication of points in the square lattice Z 2 (regarded 
as points in the complex plane) by 1 + i produces a similar sublattice of index 2. All our 
sublattices will be constructed by generalizing this remark. 

We will make use of five types of integers: Z, the ordinary rational integer s; Q, the 
ring of Gaussian integers {a + bi : a,b G Z}, where % = \J — 1; J, the ring of Eisenstein 
integers {a + buj : a,b G Z}, where u = e 2nl ^ 3 ; H , the ring of Lipschitz integral quaternions 
{a + bi + cj + dk : a,b,c,b G Z}, where k are the familiar unit quaternions; and Hi, the 
ring of Hurwitz integral quaternions {a + bi + cj + dk : a,b,c,d all in Z or all in Z + |}. 
Other rings of integers could also be used, but these suffice for the lattices considered in this 
paper. 

If A = Z, multiplication of lattice points by £ G Z gives £Z, a similar sublattice of index 

jv = iei. 

If A = Z 2 = multiplication by the Gaussian integer = a + bi E Q gives a similar 
sublattice A' = £A of index N = a 2 + b 2 . A number N is of the form a 2 + 6 2 if and only if it 
is of the form 



2ei n n 
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2% 



(21) 



=1(4) 



©=3(4) 



where the first product is over primes congruent to 1 (mod 4), the second product is over 
primes qj congruent to 3 (mod 4) and ei, fa and gj are nonnegative integers. These indices 
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7 




. 15,4 


+ 4,8 
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+ 3,7 




, 16,8 


S 8 ' 8 
8 


+ 4,4 


+ 4,0 


+ 0,3 


+ 3,0 


3 


, 3,10 




* 48 ' 4 


+ 8,0 


+ 0,4 


8°°'° 



+ 0,2 


+ 6,0 






+ 1,12 


1 


+ 1,0 


+ 0,1 


+ 2,0 


+ 2,2 


6 






+ 1 ' 5 




+ 5,0 


2 


+ 2,6 


+ 11,2 




5 

X 


+ 9,1 


5 


:^ 


+ 2,9 




6 




O 
9 




, 10,5 



10 
9 

-* 








-4 




-2 









2 





o 
11 



Figure 3: Labels generated by the algorithm. 



are the numbers 

1, 2, 4, 5, 8, 9, 10, 13, 16, 17, 18, 20, . . . (22) 



(Sequence [AT48T| of @). 



If A = A 2 = J i the planar hexagonal lattice, multiplication by the Eisenstein integer 
£ = a + bu E J gives a similar sublattice A' = £A of index N = a 2 + ab + b 2 . A number N is 
of the form a 2 + ab + b 2 if and only if primes congruent to 2 ( mod 3) appear to even powers. 
These indices are the numbers 

1, 3, 4, 7, 9, 12, 13, 16, 19, 21, 25, 27, . . . (23) 



(Sequence [A3136| of 



It is shown in [[| that the above conditions are also necessary: if Z, Z 2 or A 2 has a similar 
sublattice of index N then N must have the form described in the preceding paragraphs. 

For the lattices A = Z 4 , Z 8 , Z 12 , ... ,D± and E 8 a necessary condition for the existence of 
a geometrically similar sublattice of index is that A^ should be of the form m L l 2 for some 
integer m, where L is the dimension. This condition is also sufficient, since such sublattices 
can be obtained by writing m = a 2 + b 2 + c 2 + d 2 , regarding A as a sublattice of Hq^ 4 , and 
multiplying A on the left or on the right by the quaternion £ = a + hi + cj + dk. Left and 
right multiplications in general give different sublattices. In the case of and E s we may 
also multiply by Hurwitz integral quaternions to obtain further similar sublattices. 

Odd-dimensional lattices of dimension greater than 1 are less interesting. For a lattice A 
of odd dimension L has a geometrically similar sublattice of index A^ if and only if A^ is an 
L-th power, say m L , and sublattices of this index can be obtained by scalar multiplication 
of A by m (see H). 
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The norm of a quaternion £ = a + bi + cj + dfc is = a 2 + b 2 + c 2 + <i 2 where the bar 
denotes quaternionic conjugation. If £ belongs to one of the above rings, the index of the 
sublattice £A (or A£) in A, [A : £A], is equal to (££) L / 2 , where L is the dimension and the 
bar is complex or quaternionic conjugation as appropriate. 



B Clean sublattices 

In dimension one, the sublattice £Z is clean if and only if £ is odd. 

Reference || gives necessary and sufficient conditions for a similar sublattice of any two- 
dimensional lattice to be clean. In particular, the sublattice £Z 2 (£ = a + ib) is clean if and 
only ii N = a 2 + b 2 is odd. These indices are obtained by setting e\ = in (pl|): 

1,5,9,13,17, 25,29,37,41,45,... 

(Sequence |A57653| ). 

The sublattice (£ = a + bu) is clean if and only if a and b are relatively prime. It 
follows that A2 has a clean similar sublattice of index N if and only if N is a product of 
primes congruent to 1 (mod 6). These are the numbers 

1, 7, 13, 19, 31, 37, 43, 49, 61, 67, . . . (24) 

(Sequence |A57654| ). 

The existence of clean sublattices in dimensions greater than 2 was not considered in |J. 
We can give a fairly complete answer for the lattices Z L , L > 1. 

Theorem IV. 1 Suppose L > 1 and Z L has a geometrically similar sublattice A' of index 
N. Then A' is clean if and only if N is odd. 



Proof. (If) Let A' = </>(Z L ), where is a similarity, and let A" = 0~ 1 (Z L ). If multiplies 
lengths by c\ (as in ([I])) then N = c\ . Suppose N (and hence Ci) is odd. Let A' have 
generator matrix K, with KK tr = c\Il = ral^, where m = c\. (In the notation of ([!]), 
XJ\ = c\K\ = K.) Since A' is a sublattice of Z L , the entries of K are integers. Then A" has 
generator matrix K~ l = —K tr . 

We must show that there are no points of Z L on the boundary of the Voronoi cell of A', 
or equivalently that there are no points of A" on the boundary of the Voronoi cell of Z L . 

It is enough to consider just one face of the Voronoi cell of Z L , say that consisting of the 
points P = (|, f , f , . . . , ^), where \xi\ < 1 for 2 < i < L - 1. If P G A" then there is a 
vector u = (ui, . . . , u^) G Z L such that 

P = —uK tr . (25) 
m 

Equating the first components we get that 

11. 

- = — times a vector with integer entries . 
2 m 



Since m is odd this is impossible. 
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(Only if) Suppose N (and hence c\ and m) is even. We claim that all the vertices of the 
Voronoi cell for Z L (i.e. all the deep holes in Z 4 in the notation of ||) belong to A". In 
fact, (|25|) implies that u = PK. Let P be a vector of the form (±^, ±|, . . . , ±A and let 
K = (kij). From KK tr = K tr K = ttiIl we have Ylf=i ^% = 171 ^d, since k 2 j = /c^ (mod 2), 
Yli=i kij i s even (for all j). Hence PK has integer entries and is in Z L . ■ 

The following corollary summarizes our results about Z L for the values of L that we are 
interested in. Note that since Z L has no "handedness", there is essentially no difference 
between "similar" and "strictly similar" for this lattice. 

Corollary IV. 1 Z L has a geometrically similar sublattice of index N if and only if 

• N is an L th power, if L is odd 

• N is of the form a 2 + b 2 , if L = 2 

• N is any integer, if L = 4k, k > 1 

In each case the sublattice is clean if and only if N is odd. The same results hold if "similar" 
is replaced by "strictly similar". 

For D4 we have only a partial answer. 

Theorem IV. 2 If M is 7 or a product of primes congruent to 1 (mod 4) then has a 
geometrically strictly similar, clean sublattice of index M 2 . The values of M mentioned are 

1, 5, 7, 13, 17, 25, 29, 37, 41, 53, . . . (26) 

(7 together with Sequence \A4613j) . 



Proof. We take our standard version of the D4 lattice to have minimal norm 2 (as in 
and generator matrix 



G 



1 


1 








-1 





1 








-1 





1 





-1 





-1 



(27) 



The four rows ^1,^2,^3, ^4 of G correspond to the nodes of the Coxeter diagram for _D 4 shown 
in Fig. 4, where Vi ■ Vi — 2 (i — 1, . . . , 4), two nodes that are joined by an edge correspond to 
vectors with inner product —1, and two nodes that are not joined by an edge are orthogonal. 

We regard D4 as a subset of H = {w + xi + yj + zk : w,x,y,z G M}, the space of real 
quaternions. Our sublattices A' will be constructed by multiplying D4 either on the left or 
on the right by appropriate Hurwitzian integral quaternions. If £ = a + bi + cj + dk G Ti 
then £.D 4 has generator matrix GL^, where 



a 
-b 
-c 
-d 



b 
a 
-d 

c 



c d 

d —c 

a b 

—b a 



(28) 



15 



t'3 




• V2 



Figure 4: Coxeter diagram for any lattice that is geometrically similar to -D 4 : there are four 
generating vectors Vi,v 2 , v 3 , u 4 satisfying v 1 ■ vi = v 2 ■ v 2 = v 3 ■ v 3 = v 4 ■ w 4 , v^ ■ Vj = —\v\ ■ V\ 
if nodes Vi and Vj are joined by an edge, and ■ Vj = (i ^ j) otherwise. 



and -D 4 £ has generator matrix GR^, where 



a 


b 


c 


d 


-b 


a 


-d 


c 


— c 


d 


a 


-b 


-d 


—c 


b 


a 



(29) 



Note that 

L^Lf = R^Rf = m/ 4 , = R^ , (30) 

where m — = a 2 + b 2 + c 2 + d 2 . 

We will show that under certain conditions £.D 4 and D 4 £ are clean sublattices. We only 
give the proof for -D 4 £, the other case being completely analogous. 

The Voronoi cell for D 4 is a 24-cell, with 24 octahedral faces HQ, [|)J. A typical face (they 
are all equivalent) is that lying in the hyperplane 

1 , 
X ■ vi = -vt ■ vi , (31) 

having center 5q = \v\ and six vertices 



<*1 = 


\{**- 


f- v 3 + , 


5 2 = 


2 ( "" 3 




5s = 




f «2 + V 3 ) , 


S 4 = 




-^3), 


h = 




f- v 2 + Vi) , 


5 6 = 


2 ( "" 2 


- Vi) 



(see Fig. |^). A point X belongs to this face if and only if it satisfies (^T|) and 
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Figure 5: Labeling for center and vertices of octahedral face of Voronoi cell for D±. 



\(X - 8 ) ■ (5i -S )\ + \(X- So) ■ (S 3 -S )\ + \(X- S ) ■ (S 5 - S )\ < - Vl ■ Vl . (33) 



Let A' = D±£, where £ is a quaternion of the form 

a a . 3 . 3 , 
£ = - + -? + -? + -k , 
2 2 2 2 



(34) 



and a and 3 are odd, positive, relatively prime integers. The norm of £ is |(« 2 + (3 2 )- Then 
we claim that A' is clean. 

To show this, we begin by computing the generator matrix for A': 



G' 



GRf 






a 





P ' 


-a-P 
2 


-a+P 
2 


a-P 
2 


-a-p 
2 


a-j3 
2 


-a-P 
2 


a+P 
2 


a-P 
2 


a+/3 
2 


-a+P 
2 


-a+P 
2 


-a-P 
2 J 



(35) 



and denote its rows by v[, v' 2 , v' 3 , v' A . We will similarly use primes to denote the center (S ; ) 
and vertices (S[, . . . ,5' 6 ) of an octahedral face of the Voronoi cell of A'. From (Q) we find 
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that 

6' = ^(0,a,0,/3), 
5[ = -(a,a,P,P), 
5' 3 = ~(-(3,a,a,fl) , 

6' 6 = ~(0,a + PA-<* + P). 

We must show that it is impossible for a point X = (w, x, y, z) G D± to satisfy the primed 
versions of fl3~T|) and (fl3|), which are 

ax + f3z = ^(a 2 + p 2 ), (36) 

\aw + fiy\ + | - (3w + ay\ + \f3x - az\ < ^(a 2 + (3 2 ) . (37) 

Suppose on the contrary that (w, x, y, z) e D 4 satisfies (|36| ) and fl3~7|). From ( p6|) we have 

1 



and from (|37| 



which together imply 



z = ^-(a 2 + (3 2 - ax) (3* 



|/9x-az| < \{a 2 + /3 2 ) 



\(a-P)<x< l -(a + f3). 



So we may write rr = |(a + /x), say, where /x is an odd integer satisfying —/3 < /x < /3, and 
from (|H) 

/3 2 — a/i 

^ = 2/? ' 

which implies a/x = /3 2 (mod 2(3). Since /3 is odd, (3 2 = f3 (mod 2/5), and we conclude that 

afi = /3 (mod 2/5) . (39) 

Thus for some integer k, afi — /3 = 2k(3, and since a and (5 are relatively prime, (3 must divide 
/i. Therefore \i = ±/3. But this is impossible. For if /x = /3, x = \{a + (3), z = \{— a + /?), 
/5x — az = \{a 2 + /3 2 ), and then (|37j) implies w = y = 0, so w + x + y + z = j3^ _D 4 , since /3 
is odd. A similar argument applies if /x = — (3. 

So far we have shown that if a and f3 are odd, positive and relatively prime, then the 
sublattice _D 4 £ is clean, where £ is given by (|34|). Suppose M is a product of primes congruent 
to 1 (mod 4). From the classical theory of quadratic forms (see for example 0), we know 
that M = p 2 + q 2 with p even, q odd and gcd(p, q) — 1. We now simply set a = p + g and 

P = \p-q\- 
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It remains to discuss the case M = 7. For this we can multiply on the left or on the right 
by either of the quaternions 

+ 11.1.5, 13.3.3, 
t = — I — i H — 7H — k or — I — i H — ?H — k. 
^ 2 2 2 J 2 2 2 2 J 2 

We omit the straightforward verification that these sublattices are clean. ■ 
In the other direction we have: 

Theorem IV. 3 D± has no clean, geometrically similar sublattice of index M 2 if M is 3, 9 
or 11. 



Proof. The proof is by exhaustive search, using a computer. We produced a list of all 
vectors of norm 2M in D 4 , and from this we found all similar sublattices of index M 2 by 
finding all sets of four vectors corresponding to the Coxeter diagram of Fig. 4. Given a 
sublattice A', we compute the equations defining an octahedral face of the Voronoi cell from 
(|31|) and (|3|). Then AMPL @ and CPLEX |7| were used to verify that in every case there 
was a point of D4 on the face. ■ 

The preceding discussion has shown that the lattices Z, Z 2 , Z 4fc for k > 1 and D4 have a 
plentiful supply of clean, geometrically similar sublattices. We expect the same will be true 
of the E 8 lattice, but this question is presently under investigation. 

Finally, we remark that if A' is a clean sublattice of A and A" is a clean sublattice of A', 
then A" is a clean sublattice of A. 



C Common sublattices of A x and A 2 

We begin with a general comment. Let Ai and A2 be any two sublattices of a lattice A (they 
must have the same dimension as A but are otherwise arbitrary). Then we may form their 
intersection A n = Ai n A2 and their join A u = A\ U A2, as shown in Fig. ||. The join is 
the lattice generated by the vectors of both Ai and A 2 (and in general is not simply their 
union). From the second isomorphism theorem of group theory (e.g. f23fl ) the indices and 
determinants of these lattices are related by 

[Au : A 1 ] = [A 2 : A n ], [Ay : A 2 ] = [A, : A n ] , (40) 



detAi detA 2 = detAu det A n . (41) 

There are now in general many ways to find a "product" sublattice A s C A n with 

[A : A s ] = [A : A x ] [A : A 2 ] . (42) 

Let A be one of Z, Z 2 or A 2 , and let Ai = £iA, A 2 = £ 2 A be geometrically strictly similar 
sublattices obtained by multiplying A by elements of Z, Q or J respectively. Since these 
three rings are unique factorization rings, the notions of greatest common divisor (gcd) and 
least common multiple (1cm) are well-defined. We set £u = gcd(£i,£2), £n = l cm {£,i,£,2}, and 
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A 




Figure 6: Intersection, join and "product" sublattice of two arbitrary sublattices. 

then it is easy to see that Ay = £u A, A n = £ n A. We can also form the product sublattice 
A s = £i£ 2 A (see Fig. [?]). The indices of these lattices are given by 

[A : AJ = (£i£i) i/2 , [A : A 2 ] = (66) L/2 , 
[A : Au] = (£u£u) L/2 , [A : A n ] = (£ n £ n ) i/2 , 

[A:A S ] = [A : Ai] [A : A 2 ] . (43) 
In dimension L — 1, (^) implies that 



[A: A n ] = lcm{ [A : A 1 ], [A : A 2 ] } , (44) 



and we can take Ai cm = A n . However, if L = 2, (|44| ) does not hold in general. 

In dimensions 1 or 2, if £i and £ 2 are relatively prime (meaning gcd(^i, £ 2 ) = 1), we have 
£u = 1, £n = £i6, A = A U; A s = A n . 

Because the quaternions form a noncommutative ring their arithmetic theory is more 
complicated. For example, it is necessary to distinguish between left gcd's and right gcd's. 
Both are well-defined in Hi and also in Hq as long as at least one of the quaternions involved 



has odd norm ||, p5| . We plan to discuss this theory and its applications to the study of 
sublattices of Z 4 and elsewhere. In the present paper we will restrict our attention to 
a narrow class of sublattices, which however will be general enough to provide an adequate 
supply of sublattices for our applications. 

For Z 4 we choose two Lipschitz integral quaternions £i, £ 2 G H whose norms are odd 
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A 




A, = 66 A 



Figure 7: Join A U; intersection A n and product A s of two sublattices A 1; A 2 of A, where A 
is one of Z, Z 2 or A 2 . 

and relatively prime. For D4 we choose two Hurwitz integral quaternions 

6 = ^(l + ij + ^ + ^Hi, 

£2 = ^(l + ij + ^i + fcjel!, (45) 

where ai, a 2 , Pi, Pi are odd positive integers with gcd(ati, Pi) = gcd(a 2 , P2) = gcd((a\ + 
P 2 )/2, {al + Pl)/2) = l. 

In both cases we take A 1 = ^A, A 2 = A£ 2 and A s = A n = £iA£ 2 (see Fig. |[). Then 

[A: Ax] = (66) 2 , [A : A 2 ] = (£26) 2 , 

[A: A,] = [A : Ax] [A : A 2 ] . (46) 

For Z 4 this gives sublattices A 1; A 2 of indices M\, M 2 , where Mi and M 2 are any two 
relatively prime odd numbers (from Corollary [I For D 4 , Mi and M 2 are any two 

relatively prime numbers from (|26|) . 

V High rate asymptotics: Details 

In this section we analyze the distortion of the asymmetric multiple description lattice quan- 
tizer at high rates. 



21 




Figure 8: Ai (resp. A2) obtained by multiplying A = Z 4 or D± on the left (resp. right) by a 
quaternion £1 (resp. £2)- 

Let A be an L-dimensional lattice with geometrically strictly similar, clean sublattices 
Ax, A 2 , A n = Ai fl A 2 , A s (as in Figure ||), with indices Ni,N 2 ,N n and N s , respectively, 
where N s = NiN 2 . It is assumed that u\, the volume of a fundamental region for A, is equal 
to unity A sequence of lattices is then obtained from the base set of lattices by scaling each 
component. Let Ai(n) = nAi, A 2 (n) = nA 2 , A n (n) = nA n and A s (n) = n 2 A s . These have 
indices N n (n) = n L N n , N s (n) = n 2L N s and N s (n) = N 1 (n)N 2 (n). 

We analyze the rate-distortion performance for the set of lattices {A, Ai(n), A 2 (n),A n (n), A s (n)}. 
However, in order to keep the notation simple, we will only use the sequence index n when it 
is necessary to avoid confusion. Thus we will write A s instead of A s (n), N s instead of N s (n) 
and so on. 

Referring to ([TT|), let 



Need to 
redefine 
the Voronoi 
sets in a 
consistent 
way usinng 
notation 
that is 

^^z^Y i manageable. 



■/ s = E p ( A )E^ii A -«<( A )H 2 - ( 47 ) 



We investigate the high-rate behavior of J s and then find the approximation for di, i = 1,2. 
The latter would also allow us to predict the asymmetry in the distortion behavior of the 
quantizer. The reader is referred to Figure |9] for the analysis. Let 

J n = ^X^IK(A) - «i(A)|| 2 P(A) (48) 

and 

J„ = Efa+7 2 )l|A- 7 ' a ' (A) + w(A) || 2 f(A). (49) 

Thus, 

J s = J Sl + J S2 . (50) 
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Figure 9: Relationship of edgelength and distortion. 



Under the assumption that A n is fine enough for P(A) to be considered a constant over V\ n 
we obtain 



J si 



7i72 



71 + 72 -VeAn 



^P(A') ll«i(A)-« 2 (A)|| : 

AeV An (A') 



(51) 



By construction, the inner sum in ( |5lD does not depend on A'. Therefore, taking this out of 
the outer summation and using J2\'eAn = V-^n, we obtain 



7i72 1 
7i + 72 N n 



^2 l|ai(A) - « 2 (A)|| 



(52) 



A6V An (0) 

which can be written in terms of the edge endpoints as 

7i72 1 



J, 



7i + 72 A^ n 



E E HAx - A. 



(53) 



Xi€V An (0) A 2 eV Afl (Ai) 



Here we have used the fact that each edge labels a lattice point and therefore replaced the 
sum over the lattice points by the sum over the edges. Using the Riemann approximation 



/ | \x — Ai| \ 2 dx ~ | |A 2 — Ai| | 2 z/ 2 



(54) 



for the summation in ( j53|) we obtain 

7i72 1 



J, 



si 



E 1 



7i+72iVn A ie r n (o) z/2 



Jy s( Ai) Ik" Ail l 2 ^ 



(1+2/L) 



(1+2/L) 



(55) 



The term within the brackets is G(A S ), the normalized second moment of a Voronoi cell of 
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A s (= j2 for the square lattice), v 2 = N 2 , v s = NiN 2 and N n /[Ai : A n ] = A^. Therefore 



J Sl 



7172 [Ai : A r 
7i + 72 N n v 2 



■G(A a )u 



[ G{A S ) K 



1+2/L 



7i72 



2/L 



7i + 72 

G(A s )(iV 1 iV 2 ) 2 / L . 



7i + 72 

If all the lattices in question are scaled by /?, then 



J n =P^-G(K)(N 1 N 2 )V L . 
7i + 72 



The rate of the ith description is given by 



Therefore 



Ri = Hp) - \ iog 2 (iv,) - j io g2 (/? L ), % = l, 2. 



AT. = J_ 2 Lh(p) 2- i/? % z = l,2. 

(3 L 



(56) 
(57) 

(58) 

(59) 
(60) 



The scale factor (3 2 is related to the differential entropy of the source and the rate Rq through, 

01 = 2 2h(p) 2 -2R (61) 



Using (|60|j6T|) in ( pq) we obtain 



U S1 



7i72 
7i + 72 



Rq) 



(62) 



A bound for the term J S2 is obtained in terms of p n , the covering radius of A n , by observing 
that for every A, it is possible through a suitable A n shift, to satisfy 



A 



7i«i(A) +7 2 a 2 (A) 



7i + 72 



< Pn- 



Thus we have the inequality 



J S2 < (7i + l2)pn/3 2 



(63) 



(64) 



By comparing ( ^6] ) and ( |63D we observe that J Sl = 0(ra 4 ) whereas J S2 = 6(n 2 ). Hence J Sl 
dominates J S2 and we obtain the approximation 



J 



7i72 



■G(A s )2 2h(p) 2~ 2(Rl+jR ' 2 ~' Ro) , 



7i + 72 

where Rq determines the central distortion do and is given by do = G(A)2 2( - h ^~ Ro \ 



(65) 
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The approximations to the side distortions are obtained by using Figure |S] and the fol- 
lowing analysis. The channel 1 distortion is given by 



1 = ^- E im a )- a ip 



do- 



(66) 



AGV An (0) 



Now the central approximation in the high rate analysis of the side distortions is obtained by 
using Figure [|. The main idea being that the distance AD 2 = ||A — Ai|| 2 is well approximated 
by AB 2 = ||A — Xi\\ 2 at high rate. We will formalize this notion below. 



lA-Axll 2 = AD 2 = AC 2 
AF 2 < AC 2 < AG 2 
< CD 2 < BD 2 



CD^ 



(67) 



Now, by writing BD = e (as shown in Figure ^) and using the geometry shown we obtain 
the following inequalities (note that A = (71A1 + 72A2) /(71 + 72)), 



||Ai-A 



1 - 



< ||A- Ail I 2 < IIAi-AI 



c 




2 


e 




HA!- A||. 


+ 







Ai- A I 

Therefore we can rewrite the above as 

||A 1 -A|| 2 + ||A-A|| 2 -2||A 1 -A||||A-A|| < HA-AiH 2 ^ 
||Ai - A|| 2 + ||A- A|| 2 + 2||A 1 - A||||A- A|| + ||A- A|| 2 

As ||A — A|| 2 > and summing over A G VA n (0), we obtain 

EA £ v An( o) OlAi - AH 2 - 2||Ai - A||||A - A||] < E AG y An (o) II* - Ai|| 2 < 
Eagv^o) OlAi " All 2 + 2||A - A|| 2 + 2[|A X - A||||A - A| " 

Now, by using the fact that | |Ai — A| | < p n , we obtain 

.Easv a l|Ai-A||||A-A|f 



E 



AGV An (0) 



E 



A G y An (o) 



A : -A|| 2 
Ai-Al 



1 - 2- 



<E 



AGV An (0) 



|A- Ail P < 



1 + 2 



SA gyAn (0)HA-A|| 2 | o E, e v An (o)l|Ai-A||||A-A|| 
T,xev A (0) l|Ai-A|| 2 



Ea€V a (0) l|Ai-A|| 2 



Therefore, using these inequalities we obtain the following result. 



Lemma V.l If 71 ^ 0,72 ^ 0, liniR^oo E 
Ri — i?2 = C for some constant C . 

Proof: For our sequence of lattices 



Aevwo) 



|A - Ai| 



E 



.(6* 



(69) 



(70) 



(71) 



A€V An (0) 



|Ai — AN 2 when 
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where P is a constant that depends on A s . Hence, as = Q(n l ) we obtain lim^ 
linijj^oo — nx- = — the desired result. 
Using Lemma |V.1| we can write 



2 



E l|A-A:|| 2 « E ll*."^*^ ^ H* 1 "** 11 ' <73) 
Aev An (o) Aey An (0) Wi w A£V An (o) 

at a high enough rate. Therefore the side distortions are directly related to J s which was 
calculated earlier. Hence the side distortions are approximated by 

d - w i G(Aa)2 2fe(rt 2 -2(« 1+ fl 2 - Ro ) ( 74) 

(7i +72) 2 
2 

G(A s )2 2/l(p) 2 _2(jRl+/?2_jRo) 



(7i + 72) 



This allows us to find the approximate distortion ratio to be |^ ~ (^) 2 helping us to design 
the lattice quantizer Though this approximation is asymptotic in the rate, we observed that 
it was quite good for the quantizer design illustrated in the numerical results. 
Next we examine the case when 71 = 0, 72 7^ 0. In this case we can show that, 

ch « G(A s )2 2h{p) 2- 2Rl (75) 
d 2 « G(A a )2 2fc W2- 2 ^ 1+ ^- fl °). 

Similarly the roles are reversed when 71 = and 72 7^ 0. 

Let #0 = ^^^(l + a), then R x + R 2 - R = ^±^(1 - a). Note that a is chosen such 
that a > j|rjr^~ an d therefore Ri + R 2 — Ro < min(i?i, i? 2 )- Here we can clearly see the 
tradeoff in the central distortion do and the side distortions. 



A Minimizing average distortion 

Suppose we know that the packet loss probability on channel 1 is p\ and the packet loss 
probability on channel 2 is p 2 . Then the average distortion is given by: 

D = (1 - Pl )(l - p 2 )d + (1 - p l )p 2 d 1 + (1 - p 2 ) Pl d 2 + p lP2 E[\ |x| | 2 ] (76) 

Now using the high rate approximations developed earlier, we can find the optimal ^ needed 
for minimizing the distortion. 

Claim V.l The weights which minimize ( f76|j at high rate are given by 

71 = (1 -Pi)P2 ^ 

72 (l-P2)Pi' 
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Figure 10: Side distortions for fixed sublattices with varying 71,72- 



Proof: To optimize ([75|) we use the high rate expressions given in fl75|). Using ( 1731 ) in ([T6]) 
we obtain 



D = A + B 1 ( 



7i 



'7i + 72' 



5, 



72 



k 7i + 72' 



{7i 



where A,B 1 ,B 2 do not depend on 71,72 (they depend on R ly R 2 , Rq, (3). Without loss of 



generality, we can use 71 
substituting in (T7H) we obtain 



71 



71+72 



and 72 



^ ■ Hence defining 7 = 71 = 1 — 72 and 



D = A + B ll 2 + B 2 (l - 7 ) 2 . 



(79) 



By differentiating ( |79| ) with respect to 7 and setting it to zero we obtain the given result. 
Note that this problem is convex (just differentiate ([79]) twice and we see that it is always 
positive) and hence we have obtained the minimum with respect to 7. ■ 



VI Numerical Results 



In order to illustrate the performance of the proposed quantizer, we evaluate its rate- 
distortion performance. In order to compare its performance with that predicted by infor- 
mation theory, we assume that there is an entropy (lossless) coding of the quantizer output. 
This is done for a Gaussian source with unit variance, for which the multiple description 



rate-distortion problem was solved by Ozarow p2]| . 

The example chosen is the Z 2 lattice that we described in Section |TJ. The rates are 
chosen so that R\ — R 2 = \ log 2 ( 



JA2I- 
IAi 



In Figure 1C we illustrate the tradeoff between the two side distortions by varying 71, 72- 



In Figure |TT] we have plotted a comparison of ^ with ~. 



In Figure |12| we have plotted the side distortions and compared them with those predicted 
by information theory [[ 2|]. The key observation is that the distortion performance of the 
lattice quantizer is approximately 3dB away from that predicted by the rate-distortion bound. 
This gap is due to the shaping gain that we will pick up when we go to higher dimensions 
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Figure 11: 




Figure 12: Comparison of lattice distortion with rate-distortion bound. 

and using sublattices which have Voronoi cells which are close to spherical. The Z 2 lattice 
used in this example is more for illustrative purposes and has very little shaping gain. 

VII Discussion 

In this paper we have designed asymmetric multiple description lattice quantizers. This 
source coding scheme bridges the symmetric (balanced) multiple description quantizers and 
completely hierarchical successive refinement quantizers. Though a lattice vector quantizer 
was illustrated, this scheme could also be extended to other types of source coding schemes. 
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